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Abstract 

We construct noncommutative U(l) gauge theory on the fuzzy sphere S% as a unitary 
2N x 2N matrix model. In the quantum theory the model is equivalent to a nonabelian 
U (N) Yang-Mills theory on a 2 dimensional lattice with 2 plaquettes. This equivalence 
holds in the " fuzzy sphere" phase where we observe a 3rd order phase transition between 
weak-coupling and strong-coupling phases of the gauge theory. In the "matrix" phase we 
have a U(N) gauge theory on a single point. 

1 Introduction 

A nonperturbative regularization of noncommutative gauge theory in two dimensions is 
obtained by putting the theory on a fuzzy sphere S% pp. The lattice-like spacing parameter ( 
the inverse UV cut-off ) is proportional to 1/N where N is the size of the matrix algebra. In 
this regulator the radius of the sphere provides an IR cut-off for the theory. The limit N — >oo 
is the continuum limit [2]. 

The differential calculus on the fuzzy sphere is 3— dimensional and as a consequence a spin 
1 vector field A is intrinsically 3— dimensional. Each component A a , a = 1,2, 3, is an element 
of Mat^ with N = n(L + 1). Thus U(n) symmetry will be implemented by U(N) unitary 
transformations. On the fuzzy sphere S\ +l it is not possible to split the vector field A in a gauge- 
covariant fashion into a tangent two-dimensional gauge field and a normal scalar fluctuation. We 
can only write a gauge- covariant expression for the normal field $ as $ = ^{x a A a + A a x a + 
where x a are the coordinates on fuzzy S\ +1 defined by x a = L a /y/c2 with C2 = -§(-§ + 1). L a 
are the generators of SU (2) in the irreducible representation -|. In the continuum limit L — >oo 
the normal scalar field reduces to $ = n a A a . 
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The U(n) action on the fuzzy sphere S\ +l reads ( with the identification X a = a(L a ®l n +A a ) 
and with the normalization Tr^l^ = N ) 



S = N 



^Tr N [X a , X b f + ^ €abc TrX a X b X c 
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+ ^Tr N {Xl-a*c 2 y. (1) 



In the continuum limit this action becomes ( with the gauge coupling constant defined by 
1 

N 2 a 4 



9 2 — AT2«4 ) 



1 f dCl ^ 9 1 f dVt r „ , i r . , , , 1 2m 2 /• dfi , 9 
S =4?i iW ^'■^ b A c --[K,A i ]A c \ + - r j-tr^. (2) 

If we take the limit m — >oo ( in other words we set $ = ) then the action will reduce to a 
2— dimensional pure U (n) gauge theory 

In this last equation A a should be understood as a 2— dimensional gauge field ( in other words 
satisfying the constraint n a A a = ) and C a = —ie a b c nbd c . 

The model (pQ) with m = was obtained in string theory limit in [3]. It was shown that 
it corresponds to a gauge theory on the fuzzy sphere in [3]. It was studied numerically in [8]. 
The model with m = and m — >oo and for U(l) groups was studied in one-loop perturbation 
theory in [7j. 

The basic prediction coming from the Monte Carlo study [S] is that noncommutative U(l) 
gauge theory in 2 dimensions on the fuzzy sphere ( given by the above action ) behaves in 
the large N limit like a commutative U(N) in two dimensions on a lattice [6]. This is true at 
least in the so-called "fuzzy sphere" phase of the model for large values of the gauge coupling 
constant. Indeed we observe in the simulation that this phase splits into two distinct regions 
corresponding to the weak and strong coupling phases of the gauge field which are separated 
by a third order phase transition. This transition seems to be consistent with that of a one- 
plaquette model [10]. It seems that deep inside the "fuzzy sphere" phase the model can still be 
understood as a U(l) gauge theory on the sphere whereas in the "matrix" phase it is a U(N) 
gauge theory on a single point. 

The aim of this article is to construct models of noncommutative gauge theory on the fuzzy 
sphere which describe this third order phase transition although in general these models will 
coincide with ([1]) only in the continuum limit. An alternative approach to gauge theory on the 
fuzzy sphere is given in [9]. 



2 The Model 



The main idea is to reparametrize the gauge field on S% in terms of a single matrix W 
which contains all the tangent degrees of freedom. This W we call the " fuzzy link variable" . 
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Thus we will have the coordinate transformation (A l7 A 2 , A 3 ) — >(W,<&). Let us introduce the 
2Nx2N idempotent 

7=1(1 + 2^) , 7 2 = 1 (4) 

where a a are the usual Pauli matrices. It has eigenvalues +1 and —1 with multiplicities N + 1 
and N — 1 respectively. We introduce the covariant derivative D a = L a + A a through a gauged 
idempotent 7d as follows 

Id = 7^ , 7 = ^(1 + 2a aJ D a ) , f = 1 + + ^e abc a c F ab . (5) 

In above F ab = i[D a ,D b ) + e abc D c = i[L a ,A b ] - i[L b ,A a ] + e abc A c + i[A a ,A b }. Clearly 7 D 
has the same spectrum as 7. Thus there exists a U(2N) unitary transformation U such that 
lD = U-iU + . The crucial observation is that if we let TeU(N ~+l)xU(N -1) then UT^T + U + = 
U^fU + since T7T 4 " = 7. In other words U is an element of the d^— Grassmannian manifold 
Gm = U(2N)/(U(N + l)xU(N — 1)) and hence it ( or equivalently 7^ ) contains the correct 
number of degrees of freedom which is found in a gauge theory on the fuzzy sphere without 
normal scalar field. Indeed d N = AN 2 - (N + l) 2 - (N - l) 2 = 2N 2 - 2. Thus 7 D can be 
identified with the equivalence class [U] = {UT,TeU(N + 1) x U(N - 1)}. 

This fact can also be seen by writing down the measure in terms of U explicitly. A direct 
calculation gives Tr 2N (d'j£)) + (d'yD) = 8 Y,i=i~ Y,'j=N+2(dU + U)ij( y U + dU)ji. It is obvious that 
only d]y = 2N 2 — 2 degrees of freedom of U are involved. In the large A" limit we can see from 
(EJ) that 7£> — >7 and hence U — >l2N- Thus in this limit we can write U = 1 + O(j^) and 
as a consequence 

Tr 2N (d lD ) + (d lD ) = ^Tr N _ 1 (dA 12 ) + dA 12 + O(^). (6) 

In above A 12 is the off-diagonal upper block of A and A^ is the off-diagonal lower block of A. 
A12 is an (N + 1) x (N — 1) matrix while A+ 2 is an (N — 1) x (N + 1) matrix. We have then in 
the limit the measure 

f d lD = [ d[U]ot [ dA 12 dA+. (7) 

Integration over G^ means that we integrate over all idempotents 7^ which have A/" + 1 eigen- 
values equal +1 and N — 1 eigenvalues equal —1. A general idempotent with this property can 
be parametrized in terms of a unitary matrix U as U r jU + . Thus ( in the basis in which 7 is 
diagonal with 1jv+i in the first block and — 1n-i in the second block ) we integrate only over 
the dN = 2N 2 — 2 degrees of freedom of the unitary matrix U which correspond to jd- As we 
will see only these degrees of freedom will effectively appear in the action. 

A general idempotent with the above property can also be parametrized in terms of 3 
hermitian matrices D a as in equation (jSJ). Indeed we can check that we have ( with 7^ given 
in terms of D's by equation (jSJ) ) 

Tr 2N (d lD ) + d lD = -^Tr N (dD a ) + dD a - ^Tr N (d^) + d^ + O(^). (8) 



Then in the large N — >oo limit it is obvious that we have the measure 



f dA 1 dA 2 dA 3 = f dD 1 dD 2 dD 3 oc f dj D f d®. (9) 

In the limit m — >oo we have $ = and hence we obtain instead 

f dD 1 dD 2 dD 3 5($) oc f dj D . (10) 

Thus for consistency we must set $ = in 7 2 . We can then compute ( with the covariant 
coordinates x% defined by = j?D a and eaF = e abc a c F ab ) the following expansion 

7^ = Jj + ™° - ^(^)(«xF) - ±- 3 (eaF) + ^(^)( e( xF) 2 + ^(eaF) 2 + 0(^.(11) 

Let [/ and be 2iV x 2N unitary matrices which are in the Grassmannian manifold and 
let us consider the following path integral 

Z[j,X]= / d[U] / d[V] exp(-Sp) , S P = — + -Tr 2N U 1 U + V 1 V + + h.c. (12) 

JGn JGn A A 

As we have explained integration over Gn means that we integrate only over the 2d^ = 2(2N 2 — 
2) degrees of freedom of the unitary matrices U and V which will effectively appear in the action. 

The matrix U r yU + can always be parametrized as U r yU + = 7^ where 7^ is the idempotent 
given by ([3}. For small U near the identity given by U = 1 + i4 + ... we can show that 7/3 
( and hence the action ) will only depend on the off-diagonal blocks A12 and Af 2 of A. We 
will only integrate over these components in the measure. Similarly the matrix V^V + can be 
parametrized as V r yV + = 75 where 7^ is another idempotent given by equation (jSJ) but in 
terms of a different covariant derivative B a = L a + A' a with a new gauge field A' a and normal 
component $ . As before we write V = 1 + ijr + ... in the large N limit and again the action 
will only depend on the off-diagonal blocks A' 12 and (A' 12 ) + of A'. We will only integrate over 
these components in the measure. 

The action Sp is invariant under all U(2N) unitary transformations of the form U — >gU, 
V — >gV which is only possible beacuse of the doubling of gauge fields. The idempotents 
7d and 7^ will therefore transform as 7£> — >g , ~fDg + and 7b — >giB9 + respectively. Thus the 
measures djo and djp are also U(2 N)— symmetric and therefore the full path integral ffl2l) is 
also U(2 N) — symmetric. 

Next we compute the action Sp in the configuration W = JdJb- Let G ab and xf be the 
curvature and the covariant coordinates associated with the covariant derivative B a . We obtain 
after a short calculation 



Sp = ik TrN 



- 32c 2 + 2N\x» + x B a f + j^[D a - B a , D b - B b ]{F ab + G ab ) - — {F ab + G ab ) 



+ + G 2 J - ^x°F cd x* G cd - ± x *F cd x°G cd 

+ l^ x °(eF) a x*{eG) b + jL x °(eF) b x*(eG) a + l^x B a {eF) a x° (eG) b + ^(eF^x^eG), 



/ r / ) /? ; r /? Dl\ / / / / ~s / «\ \ I 
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V o , rf,*f] + ( (eF).(eF) 6 + (eG) a (eG) b j + 0(— ) 



(13) 



These are the first few terms of the U(l) action in this fuzzy plaquette model. The plaquette 
variable is the unitary 2N x 2N matrix given by W = JdJb- This fuzzy plaquette U(l) model 
is more complicated than the original fuzzy U(l) model (pQ) although the two actions (jTJ and 
(j!3p have the same continuum limit as we will now show. 

We show this result in two steps. First we compute the continuum limit of the above classical 
action then we integrate out one of the gauge fields. We will see in particular that the effect 
of the first non-trivial term 2N 2 Tr^{x^ + x^ ) 2 is such that the path integration over B a is 
dominated by the configuration B a = —D a . Indeed in the continuum large N limit in which 
we keep N 2 X fixed we can use in the classical action Sp the limits x^,x^ — m a where n a are 
the global coordinates on the sphere and hence Sp reduces to 

F a b + G ab — QF ab G a b + 2(e a b c n c F afe )( 

^abcTl cGab) 

In above we have used the fact that on the sphere {e a b c n c F a b)(e a b c n c G a b) = 2F a bG a b and * stands 
for the star product on the fuzzy sphere which still appear in the second term. The path integral 
over the B a is clearly seen to be dominated by the configuration B a = —D a l . Therefore we 
obtain the action ( modulo a constant term ) 



Si 



62c 2 21\' 



A 



A 



— ( X D + x B\i + _L_ 
4tt 1 a + %ah + N 2 X 



47T 



Sp = ^Tr 2N W + b ,= ^-J^ b . (15) 

This action is essentially the U(l) action obtained from (jTJ) ( see equation (jSJ) ) provided we 
make the following identification 

16 _ 1 _ a 4 _ a 4 
N 2 X ~ Ag 2 ~ 4 ~ AN 2 ' [ } 

Hence the fuzzy U(l) action ([1]) with fixed coupling constant a corresponds in this particular 
limit to the fuzzy £7(1) plaquette action ( TTBl) in the weak regime A — >0 and agreement between 
the two is expected only for weak couplings ( small values of A or equivalently large values of 
at ). This is what we observe numerically [5]. 



3 The solution 

The path integral (ff2l) is invariant under all U(2N) unitary transformations of the form 
U — >gU, V — >gV or equivalently 7d — > 5'7d5 ,+ , 1b — >91b9 + - This symmetry can be fixed as 
follows. We perform the following U(2N) transformation U — >U = gll with g = V + . For 
fixed V we obtain the measure J Gn dll = J Gn dll' and the action Sp will only depend on U . 

^^This means in particular that the unitary matrix W = jdJb approaches in the quantum theory — I2.1v in 
contrast with the classical limit W — >l2JV- 
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The integral over V can be done and one ends up with the path integral ( by dropping also the 
primes ) 



r /AN 2 N 

Z[j, A] = J g d[U] exp - -Tr 2N ( 1 U 1 U + + h.c 

(AN 2 N \ 

d^ D exp — r Tr 2A r(77 D + h.c) . (17) 

G N V A A / 

This is still invariant under U(N + 1) x U(N — 1) unitary transformations U — >gll. This is 
the gauge symmetry we want on the fuzzy sphere with N x N matrices. Let us recall that we 
want a U (N) gauge symmtery on the fuzzy sphere but in this construction where we are using 
2N x 2N matrices it is only natural to obtain a gauge group which is twice as large. Remark 
also that we do not get precisely U(N) but the groups U(N + 1) and U(N — 1). In the large 
N limit this becomes unimportant. 

Apart from the restricted integration which is performed on the Grassmannian manifold Gn 
the path integral (ITT)) looks very much like a U(2N) one-plaquette model. Indeed the matrix 
W = 77d is a unitray 2N x 2N matrix. Let us then write (ITTj) in the equivalent form 

Z[ 7 , A] = jdW J d lD 5(W - 77^) exp - jTr 2N (W + h.c)) 

r /AN 2 N \ 
= J dW I(W) exp f — —Tr 2N (W + h.c) J . (18) 

We can therefore approximate the partition function of noncommutative U(l) gauge theory on 
the fuzzy sphere given by (1121) with the following path integral 

r /AN 2 N \ 

Z[ 7) A, M 2 ] = J dW I(W, M 2 ) exp (—- - —Tr 2N {W + h.c) J 

I(W,M 2 )=J g d 7D exp(-M 2 Tr 2N (W - llD ) + (W - llD )y (19) 

Instead of the delta function in (fl8|) which implements the constraint W—'yyr, = we add a new 
term to the action with large positive coupling constant M 2 which implements this constraint 
only approximately. As it turns out the integral over Gn in (|T9l) is much easier to compute ( 
at least in large N ) than the original corresponding path integral I{W) in (j!8j) . 

We remark that by dropping the requirement that M 2 must be large we get a model with an 
enlarged phase space. Indeed the space of parameters of the theory becomes the two dimensional 
quadrant A>0, M 2 >0 instead of the original positive real line. An alternative way of thinking 
about (JT9~j) is as follows. The one-parameter family of models Zm = Z[y, A, M 2 ] interpolate 
between the model Z which corresponds to a true U (2N) one-plaquette model and Z^ which 
is the path integral of the noncommutative U(l) gauge theory on the fuzzy sphere given by (|18|) 
or equivalently (fl2!) . In other words holding N fixed and taking the limit M 2 — >oo reproduces 
the original path integral (Tl8l) . 

It is hopped that the models Zm with intermediate large values of M 2 will capture some 
of the most important features of noncommutative U(l) gauge theory on the fuzzy sphere ( at 
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least in the "fuzzy sphere" phase ) seen in the Monte Carlo study of the action flTJ). In particular 
we will consider in the following the double scaling limits M 2 , N — s>oo keeping M 2 /N fixed. By 
construction the models Zm coincide in these double scaling limits with the limit iV — >oo of 
( I18p . So these models have the correct continuum limit and hence they are also path integrals 
of noncommutative U(l) gauge theories on the fuzzy sphere. However quantum mechanically 
these models for finite fixed values of the ratio M 2 /N are found to behave differently from 
( TTgj) which corresponds to the ratio M 2 /N — >oo. As we will see the limit M 2 /N — >0 is also 
different. 

The matrix 77^ is a unitray matrix which in the continuum large N limit tends ( by equation 
([5]) ) to the identity matrix 1 2 jv- Indeed in this limit both 7 and 7^ approach the usual chirality 
operator 7 = n a a a and hence 77,0 — ►I2.1v- Equivalently the unitary matrix U will approach in 
the continuum limit the identity 1 2 n as U = 1 + ijj + 0(jj%). In this limit it is also seen by 
equation ([7j) that the measure J Gn d^r, = J Gn d[U] depends only on the off-diagonal blocks A 12 
and A 2 i with A 21 = A i2 . Let W\ and W2 be the diagonal blocks of W in the basis where 7 is 
diagonal. We denote by W12 and W21 the off-diagonal blocks. Thus in the large N limit the 
second line of ( 1T51) takes the form 

I(W,M 2 ) = f d-i D e- M2Tr2N ^ 2 - w+llD - WlD ^ 

JGn 

oc J d^2dK2ie- M2Tr2N ^ w - i y^ w ~ i y^ TrN + 1 ^- Ww ^ 

(20) 

There are two cases to consider. 

Case I : This corresponds to the special case M 2 = 0. By the requirement of continuity this 
point should be reached in the double scaling limits M 2 ,N — >oo keeping ^ fixed such that 
^ >0. The above integral (!20j) becomes 

I(W,M 2 ) = e -M 2 Tr 2N {w-l) + {w-l)_ (21) 

Thus we end up with the partition function 

r /AN 2 N Nil M 2 

Z[ 7 ,A,M 2 ](x/rfH/ exp ( Tr 2N (W + h.c)) , — = --— (22) 

J V Am Am / Am A iv 

Case II : This is the generic case of the double scaling limits M 2 , iV — >oo keeping ^ finite 

and fixed. Extrapolation of the results to the case ^ >oo will by definition correspond to 

the model (Tl8l) . However extrapolation of the results to M*- — >0 will not reproduce case I as 
we will see. In these double scaling limits the integral (TSUI) becomes 

I(W,M 2 ) oc e-^'^ 2 -^-^)-^^^ (23) 

In above Wi is an (N + 1) x (JV + 1) matrix, W 2 is an (N — 1) x (N — 1) matrix and W 12 , A 12 
are (N + 1) x (N — 1) matrices whereas W21, A 2i are (N — 1) x (N + 1) matrices. 
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Since in these limits M 2 is large proportional to iV we also see that we have W — >1 2 n or 
equivalently W\ — >l/v+i an d W 2 — >-1jv-i and hence W\ 2 + W 2 \ — >0 which follows from the 
identity WfW 12 + W 2 \W 2 = 0. Together with the delta function 8{W 2 \ - W 12 )8{Wf 2 - W 21 ) 
appearing in the last line of equation fl23l) we can conclude that the off-diagonal blocks W\ 2 and 
W 2 i of W will approach zero in this double scaling limit. By neglecting also the edge effects in 
the large N limit we can take W\ and W 2 to be N x iV matrices. The partition function given 
by the first line of (119]) becomes 

r /2N 2 N \ 

Z[ 7 ,A,M 2 ]oc[Z(A M )] 2 , Z(X M )= dW 1 exp - Tr>(Wi + h.c) V (24) 

J \ am am / 

Therefore in these double scaling continuum limits in which M 2 , N — >oo keeping M 2 /N fixed 
we can set I(W) in the first line of equation (|T9l) equal to 1 and replace W with the matrix 
obtained by taking the diagonal parts W\ and W 2 to be two arbitrary independent unitary 
N x N matrices while allowing the off-diagonal parts Wi 2 and W 2 i to go to zero. In this 
approximation we can drop the restriction that these unitary matrices W\ and W 2 must be 
close to Ijv since we have replaced the coupling constant A with Am- 

The path integral of a 2— dimensional U(N) gauge theory in the axial gauge A\ = on a 
lattice with volume V and lattice spacing a is given by Z(\ M ) v/a2 . Therefore we can see that 
the partition function Z(7,A,M 2 ) of noncommutative U(l) gauge theory on the fuzzy sphere 
is proportional to the partition function of a 2— dimensional U(N) gauge theory in the axial 
gauge A\ = on a lattice with two plaquettes. This is true at least in the weak coupling region 
of the phase space. Indeed since Am>0 we must have 

, N , s 

(25) 

The matrix phase The fact that the one-plaquette model ( J24l) can be defined for all values 
of the coupling constant Am in the range [0, oof leads to the above restriction on the allowed 
values of the gauge coupling constant A. This can be understood as follows. The partition 
function Z^j, A, M 2 ) corresponds to noncommutative U(l) gauge theory on the fuzzy sphere 
only for values of the coupling constant A below the upper value -™ where we can make sense 
of this model as a one-plaquette model. For A>-^j the model Z{pj, A, M 2 ) is presumably in the 
so-called "matrix" phase where there is no an underlying stable sphere as a spacetime and the 
gauge theory is defined on a single point. 

This picture is consistent with our previous results OCT]. Indeed we found in the one- loop 
calculation as well as in numerical simulation that the model ([I]) undergoes a first order phase 
transition from the "fuzzy sphere" phase to a "matrix" phase where the fuzzy sphere vacuum 
collapses under quantum fluctuation. The fuzzy sphere phase is defined only for values of the 
coupling constant A such that ( m 2 is the mass parameter appearing in (JTJ ) 

A<Af^°> = ^(m 2 + v / 2-l). (26) 
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The one-plaquette 3rd order phase transition In this section we solve the models given 
by (j22p and (I24p in the large N limit. The only difference between the two models lies in the 
fact that (l2~2~j) is a U(2N) plaquette model while (T2"l|) is a U(N) plaquette model and hence we 
will concentrate on the detail of (liHj) which describes noncommutative U(l) gauge theory on 
the fuzzy sphere. 

We can immediately diagonalize the matrix W\ in (T23]) by writing W\ = TDT + where T is 
some U(N) matrix and D is diagonal with elements equal to the eigenvalues 9i + tt of W± 2 . In 
other words = 8ij(9i + tt). The integration over T can be done trivially and one ends up 
with the path integral 

N 



= / H*>« ^ ■ ^ = ^ + T 2 -Ecos^ + -^Eln(sm^)'. (27) 
•/ Am i= i 2iV ^ V 2 / 



The last term is due to the Vandermonde determinant. In the large N limit we can resort to 
the method of steepest descent to evaluate the path integral Z(X) . The partition function will 
be dominated by the solution of the equation = which is a minimum of the action SV 
In the large N limit we introduce a density of eigenvalues p(8) which is positive definite and 
normalized to one. Thus sums will be replaced by J2i — N J d9p{9). The saddle point solution 
must satisfy the equation of motion 

sin 9 = [ rfrp(r) cot ^JlL, (28) 



A 



M 



By using the expansion cot = 2 Y^=i ( s i n n @ cos nT ~ cos n & nr ) we can solve this equa- 
tion quite easily in the strong-coupling phase and one finds the solution [in] 



p(0) = -L + -L C os0. (29) 

Z7T TCAm 

However it is obvious that this solution makes sense only where the density of eigenvalues is 
positive definite. This density of eigenvalues p is positive definite for all values of 9 iff Aj\/>2. 
The critical value of the coupling constant Am is then seen to occur at Am* = 2. In terms of 
the coupling constant A we obtain the critical value 

K = ^w- (30) 



1 -I- 2- 

N 



From ([25]) we then see that for the weak coupling phase we have \<\*<jj?- In the strong 
coupling phase A"<A<-ga. 

In the weak regime we obtain the solution pU] 



^ ) = ^ COS 2V~ Sm 2' Sm Y = V^- (31) 

2 The addition of the angle 7r reflects the fact that the unitary matrix W approaches in the quantum theory 
— I2N and not I2N for very weak couplings. See the footnote on page 5. 
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Again the critical value of the coupling constant Am is seen to occur at Am* = 2 when the 
critical angle approaches ir. At this point the eigenvalues 9 of W\ will fill the whole unit 
circle. 

For very strong couplings the density of eigenvalues (129]) becomes a uniform distribution 
while in the very weak coupling region ( which corresponds to very small angles 6 ) the density 
of eigenvalues f[3"Tj) reduces to Wigner semicircle law, viz 

1 



P (6) = —p\ M -6 2 . (32) 

We can use this last equation to compute the free energy and specific heat for small values 
of the coupling constant A with excellent agreement with the simulation results of (JT]). The free 
energy is given by 

F„ „ S N 8 2 r e * „ fO* r e * , , _ (9-a^ 2 



tJ = 2x| = T^-f- /* d9p(9)9 2 + r d9p{9) r dap{a)\n 
N 2 JM A A / Am j-e* J-e, J~e, 

= ^-1 + ln^ + c, (33) 
Am * 



Sn 
N 

one-plaquette models contributing to F in accordance with ff24|) . In order to compute the 

F 
' T 



The factor of 2 which multiplies comes from the fact that we have two identical independent 

to F in accordance with 
specific heat we implement the scaling transformations F — >^ and Am — >T\m- The specific 



heat is then defined by 



an^L- (34) 



A straightforward calculation yields the very simple result 

w = i- < 35 > 

This is precisely what we see numerically in the fuzzy sphere phase. The value 1 emerges 
from the fact that we have two plaquettes. As it turns out this result is valid throughout the 
weak-coupling phase, i.e for all Am<2. 

In the regime of strong couplings the free energy and specific heat are computed using the 
distribution of eigenvalues (129|) . We find 



N 2 



2x^ = ^ + ^/ d9p{9)cos9+ d6p{6) dap{a) In sin — — 



\M am J-n 

4 4 i 



4 4 /" 7r 

+ —5 — I — 7T-7T- / dx(ir — x) cos 2x ln(sinx) 2 + C2 
Am Ki ttAm-Zo 
4 2 

= ^~ + TT + C2- (36) 
am am 

The specific heat in this phase is therefore given by 

(37) 



C v , 2 2 



N 2 V A 



M 
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The phase diagram From equation (!30|) we can immediately see that for every value of the 
parameter M 2 /N the model (1241) undergoes a third order phase transition from strong-coupling 
phase to weak-coupling phase consistent with the ordinary one-plaquette third order phase 
transition observed in 2 dimensional gauge theory on the lattice [10J. However in this case this 
transition occurs at the value A" of the gauge coupling constant A which becomes smaller as 
we increase M 2 /N until it vanishes for M 2 /N — >oo. By definition the double scaling limit 
in which the fixed ratio M 2 /N is taken to infinity corresponds to the model (fTSj) . Indeed the 
partition function (fT8l) is obtained from the partition function Z[j, A, M 2 ] given by (fl~9l) in the 
limit in which we take M 2 — >oo first and then N — >oc. Hence we can immediately conclude 
that the model (fIBl) does not undergo the above one-plaquette phase transition since A" — >0 
when M 2 /N — >oo and furthermore this model according to the restriction ( 1251) exists mostly 
in the matrix phase. 

Thus the model (124"!) with finite values of the ratio M 2 /N behaves as an ordinary 2d lattice 
gauge model with a reduced one-plaquette critical value whereas in the limit M 2 /N — >oo the 
3rd order one-plaquette phase transition is completely removed from the model. This limit can 
be thought of as a way of regularizing the one-plaquette transition in the lattice model. 

However for M 2 /N — >0 the above computed critical value A" = 2 is not the correct critical 
value at which the third order phase transition should occur. The model (fT9|) with the fixed 
ratio M 2 /N such that M 2 /N — >0 is a U{2N) one-plaquette model which should be described 
by the partition function (|22|) and not by (I24j) . By going through the same steps which led to 
( 13H|) we obtain for the model (1221) the critical value Am* = 1 or equivalently 

1 M 2 M 4 

* = TP£ = 1 -7T + °(^- (38) 

In above we have used the fact that the partition function (1221) should really be used only for 
very small M 2 /N. Clearly for M 2 /N — >1 the critical value A' — >0 and then it becomes negative 
for M 2 /N>1 and hence the above formula should only be valid in the range M 2 /iVe[0, 1[. In 
other words the 3rd order one-plaquette phase transition occurs at the values X[ of A for all 
M 2 /N&[0, 1[ whereas for M 2 /N&[1, oo[ the transition should occur at the values A". The value 
M 2 /N = 1 seems to demarcate a discontinuity in the phase diagram. Below M 2 /N = 1 the 
model (fl~9l) describes an ordinary U(2N) gauge theory on one plaquette whereas above this value 
the model becomes the U (N) gauge theory on two plaquettes given by (12~4"1) which corresponds 
to noncommutative gauge theory on the fuzzy sphere. The numbers of degrees of freedom of 
the model below and above this critical mass M 2 /N = 1 are different. Above M 2 /N = 1 we 
have 2N 2 degrees of freedom whereas below this value we have 4N 2 degrees of freedom. 

Thus the model fflHj) with M 2 /N = 1 is seen to be the first noncommutative U(l) gauge 
model on the fuzzy sphere Sfj which we will obtain as we increase M 2 /N from to oo. Indeed 
this model has the correct number of degrees of freedom given by 2N 2 and it can be reduced 
to (124"1) . The model with M 2 /N = 1 is also the noncommutative model with the largest "fuzzy 
sphere" phase in the sense of (I25I) hence the effects of the "matrix" phase are the weakest in this 
case. As a consequence the estimation of the critical point of the 3rd order phase transition is 
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The 3rd order one-plaquette phase transition 



2.5 



equation (25) 
equation (30) 
equation (38) 



matrix phase 




Figure 1: The weak-coupling and strong-coupling phases of the noncommutative U(l) gauge 
theory on the fuzzy sphere are the regions of the weak-coupling and strong-coupling phases 
with M 2 /N>1. The regions with M 2 /N<1 correspond to an ordinary U(2N) gauge theory. 
The number of degrees of freedom of the model above the critical mass M 2 /N = 1 is 2N 2 
whereas below M 2 /N = 1 we have 4N 2 degrees of freedom. The fuzzy sphere phase consists 
of the weak-coupling and strong-coupling phases with M 2 /N>1. In this phase the theory is a 
U (N) gauge theory on two plaquettes. 



most reliable using the model (fl~9l) with M 2 /N = 1. At M 2 /N = 1 the critical value A" is given 
by ^* = 2/3 = 0.66. This computed critical value leads to the critical value of the coupling 
constant a ( from equation (TIB"]) ) 

64 

at = — = 96 a* = 3.13. (39) 

K 

This is to be compared with the observed value a* = 3.35±0.25 seen in Monte Carlo simulation 
of the model ([I]) with the Metropolis algorithm. The agreement is very good. Thus it seems 
that the two noncommutative U(l) models given by (CQ) ( with large m 2 ) and ( [191 ( with 
M 2 /N = 1 ) are quantum mechanically equivalent in the fuzzy Sfj phase. 



4 Conclusion 

We constructed in this article a one-parameter family of noncommutative U(l) gauge theory 
partition functions Z M = Z(j, \,M 2 ) on the fuzzy sphere Sff given by equation (fl9|) . The 
idempotent 7 is the chirality operator and A is the gauge coupling constant. In the limit 
M 2 — >0 the model becomes a U(2N) one-plaquette model whereas in the limit M 2 — >oo we 
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get the generalized U(N) one-plaquette model defined by (fTBl) . In the double scaling limits 
M 2 , N — >oo keeping M 2 /N fixed >1 we can show that quantum noncommutative U(l) gauge 
theory on the fuzzy sphere Sfj is equivalent to a U(N) nonabelian Yang-Mills theory on a 
two dimensional lattice with two plaquettes according to ([21]). Thus the model ffT9l ( with 
M 2 /N = 1 ) will undergo in the fuzzy sphere phase a 3rd order one-plaquette large N phase 
transition between weak-coupling and strong-coupling phases around the critical value (13"9l 
which is consistent with the observed 3rd order transition point of ([I]). In the "matrix" phase 
we will have a U(N) gauge theory on a single point. 

Another evidence for this equivalence comes from the calculation of the specific heat which 
seems to agree with the simulation results of (JTJ) both in the weak-coupling and strong-coupling 
phases up to the sphere-to-matrix first order transition where the whole spactime ( the sphere 
) collapses. In particular the specific heat was found to be equal to 1 in the fuzzy sphere weak- 
coupling phase of the gauge field which agrees with the observed value 1. The value 1 comes 
precisely because we have two plaquettes which approximate the noncommutative U(l) gauge 
field on the fuzzy sphere. In the strong-coupling region deviations between ( l37l) and the data 
coming from the simulation of (pQ) become significant only near the sphere-to-matrix transition 
point [51 E] . 
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